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Passive states are defined as those states that do not allow for work extraction in a cyclic (unitary) 
process. Within the set of passive states, thermal states are the most stable ones: they maximize 
the entropy for a given energy, and similarly they minimize the energy for a given entropy. Here 
we find the passive states lying in the other extreme, i.e., those that maximize the energy for a 
given entropy, which we show also minimize the entropy when the energy is fixed. These extremal 
properties make these states useful to obtain fundamental bounds for the thermodynamics of finite¬ 
dimensional quantum systems, which we show in several scenarios. 


I. INTRODUCTION 

Passive states were introduced in [1] as the ones obey¬ 
ing the second law of thermodynamics in the Kelvin- 
Planck formulation [2, 3]. Namely, states that can yield 
no work in a Hamiltonian process at the end of which the 
system returns to its initial Hamiltonian, H. Any such 
process can be described by a unitary operation U, and if 
we define the maximal extractable work from the system 
as 


Wniax(p) = maxtr [H {p - UpU^)] , (I) 

then passive will be the states for which Wmax = 0. The 
quantity Wmax was given the name “ergotropy” [7] . 

Although the second law is formulated for thermal 
states [4], passive states constitute a much wider class 
[5]. In fact, they consist of all states that commute with 
the system Hamiltonian and have no population inver¬ 
sions [1, 6, 7]. Thermal states enter the picture in two 
ways. First, they are the ones that have the minimal 
energy for a given entropy. Second, thermal states are 
the only completely passive states. Complete passivity is 
another fundamental notion introduced in [1], and desig¬ 
nates those states p for which p®" are passive for all n. 
The fact that only thermal states are completely passive 
is very well illustrated by the elegant result in [8] , stating 
that the asymptotically activatable work contained in a 

passive state Wact = lim ^ jg given by 

n—>00 ” 

Wact = tr(i7crp) - tr(iJTp), (2) 

where rp is the thermal state at the inverse tempera¬ 
ture /3 [4], and (3 is uniquely determined by requiring 
S{ap) = S{Tp). Here S{p) = — tr(plnp) is the quantum 
von Neumann entropy [8]. Summarizing, provided it is 
not thermal, a passive state can be activated by jointly 
processing several copies of it, and the work that can be 
extracted in the limit of infinite copies is given by (2). 

Being motivated by the above results, the main goal of 
this work is to identify and study the other extreme of 


passive states (with respect to completely passive, ther¬ 
mal states), i.e., the ones that have maximal energy for a 
given entropy [9] . We refer to such states as the most en¬ 
ergetic passive states (MBPS). We then show that, due to 
their extremal properties, these states provide useful in¬ 
formation about fundamental thermodynamic processes. 
First, from their definition, it naturally follows that the 
MBPS have maximal activatable work content, see (2). 
Another motivation for our study is that while thermal 
states, when used instead of UpU^ in (1), provide an up¬ 
per bound on the extractable work, the MBPS provide a 
lower bound. From a methodological point of view, this 
gives a practical tool to estimate the usefulness of a given 
state from the perspective of average work extraction. 

Akin to thermal states, the MBPS have a rather gen¬ 
eral characterization and are also monotonic with respect 
to entropy. They constitute a one parameter family, and 
take a particularly simple form, 


P = 


z=0 


1-A 




2 = 0 


( 3 ) 


where \ei) are the energy eigenvectors, and e^+i > e^. 
That is, the state is (at most) a mixture of two projec¬ 
tors onto subspaces of states with energies lower than a 
given value. These states are known as 0-canonical states 
[10], and are exactly the passive states related to micro- 
canonical states. This gives a new meaning to this rarely 
used concept. 

The MBPS also allow us to quantify how energetically 
different passive and thermal states can be. Quite re¬ 
markably, it turns out that although the MBPS can de¬ 
viate significantly from thermal states for different spec¬ 
tra, we give evidence that the MBPS of, e.g., many-body 
systems with short-range interactions, behave almost as 
thermal states and have little potential for locked (i.e., 
potentially activatable) work. This makes another case 
for the universality of the thermodynamic formalism in 
the macroscopic world [3, 11-14]. 
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II. PASSIVE STATES 

Consider a process where the system remains ther¬ 
mally isolated: it can evolve according to its own Hamil¬ 
tonian H and due to external (time-dependent) helds, 
V{t). Furthermore, the process is cyclic, i.e., the ex¬ 
ternal fields are turned on and off at the beginning and 
at the end, H(0) = V(t) = 0. The corresponding evo¬ 
lution can be described by a unitary evolution U, with 
U{t) = ex^ {—i JJ" dt [H -|- V(t)]}. Since the system re¬ 
mains thermally isolated, work is given by the change of 
its average energy, W = tv {pH) — tv{UpWH), where 
H = ei|ei)(ei| (e^+i > e^) is the internal Hamiltonian 
of the system. 

In this process, work is provided (or extracted) by the 
external time dependent fields V{t). By appropriately 
choosing V {t) , we can generate every unitary operation 
H, and thus the operations considered in this context are 
equivalently all unitary operations. It follows that the 
maximal work that can be extracted from p is given by 
(1). This expression is maximized for UpU^ = ap, with 

ap = '^Pi\ei){ei\, with p^+i < Pi, (4) 

i 

where pi are the eigenvalues of p [1, 6]. In other words, 
given a state p, the (maximal) extractable work reads 
Wmax = tv{pH) - tr(crpH) [7]. 


III. MAIN RESULT 



FIG. 1. (a) The set of all passive states S and the intersection 
with the constant energy hyperplane S£ for a four dimen¬ 
sional system, (b) Entropy versus energy given an equally 
spaced Hamiltonian of 4 levels, i.e., H = diag{0,1, 2,3}. The 
shaded area corresponds to the simplex 5. The two boundary 
curves correspond to the set of thermal states (upper) and the 
most energetic passive states (lower). 


the polytope formed by this intersection, i.e., SS = {ap : 
ap G S,tv{Ha) = E}. 

The point then is to minimize the entropy function 
S{a) = — tr((T In cr) over SS. These considerations are 
illustrated in Fig. 1. Since S£ is a polytope and the en¬ 
tropy is a concave function, the minimum is achieved at 
a vertex [16], each of which have a simple form [17]. The 
vertices occur at the intersections of the energy hyper¬ 
plane with the edges of S and therefore have the form 
ap{k,l) = Xujk -b (1 — X)uji, where A = X{k,l) is deter¬ 
mined from the energy condition: 


X{k,l) 


tv{HuJi) — E 
tv{Hu:i) - tv{Hu}k)' 


( 7 ) 


In this section, for a given Hamiltonian H and entropy 
S [15], we find the passive state that maximizes the en¬ 
ergy, which we denote by a*. This maximization yields 
an upper bound on FFact in (2), 

Wact < tv{{a; - Tp,)H) = A^ax(5, E) (5) 

where S{p) = S{Tp') = S, and E = tr(iJp). It is conve¬ 
nient to first consider the complementary optimization, 
i.e., to find the passive state that minimizes the entropy 
for a fixed energy E. We will show that both optimiza¬ 
tions provide the same state. 

It is useful to introduce the following set of d linearly 
independent states: 

1 ^ 

Wfc = - ^ |ej)(ei|, l<k<d. (6) 

Any passive state can be written as a convex combination 
of such states, ap = with Qi > 0 and qi = 

1. Therefore, the set of passive states defines a convex 
polytope (in fact, a simplex) which we denote by S, whose 
vertices are given by Wk in (6). 

Within S, we are interested in the subset of states with 
constant energy, tv{pH) = E. Since the energy tv{pH) 
is a linear function, the condition tv{pH) = E dehnes an 
hyperplane, which intersects with S. We denote by S£ 


Note that for consistency tv{HuJk) < E < tr{HuJi) must 
be satisfied, i.e., the vertices must be separated by the 
energy hyperplane. In general, the set of feasible index 
pairs 2{£) = {{k,l)\tv{HoJk) < E < tv{HuJi)} will de¬ 
pend on the spectrum of the Hamiltonian and the aver¬ 
age energy E. It is however efficient to enumerate, with 
a system of dimension d requiring only to check 0{d^) 
pairs. The last step of the optimization is to minimize 
the entropy over all feasible pairs 

crp{E) = i^nap{k,l), ( 8 ) 

which can again be carried out efhciently for finite di¬ 
mensional systems. 

We denote the entropy as S*{E) = S{ap{E)). If it 
is a monotonically increasing function of A, then ap{E) 
is also a solution of the complementary optimization, 
namely maximizing the energy when the entropy is fixed. 
In the following we show that this is the case by reduc- 
tio ad absurdum. We dehne the polytope of all pas¬ 
sive states with an energy greater than or equal to E, 
S£+{E) = {ap : ap G S and tr(iJcr) > E}, whose ver¬ 
tices are those of S£ plus those of S whose energies are 
at least E. Again, the minimum of S{a) over 5£+(A), 
S{_{E), is achieved at one of the vertices. Assume that it 
is one of the with tv{HuJk) > E. Consider the passive 
state awi -I- (1 — a)ujk, with A given by A(fc, 1) in (7), so 
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that its energy is equal to E. A direct calculation shows 
that S{auii + (1 — a)uJk) < S{u!k), which contradicts our 
previous assumption. This implies that the minimum of 
S{a) over S£+ is attained on S£, which, along with the 
observation that S£+{E') C S£+{E) if E' > E, shows 
that the entropy is a nondecreasing function of E. 

In conclusion, the passive states that maximize the en¬ 
ergy for a fixed entropy, and at the same time minimize 
the entropy for a given energy, are the one parameter 
family defined by (8). This family lies on the boundary 
of the set of passive states (see Fig. 1), which are con¬ 
vex combinations of states given by (6). This suggests a 
beautiful relation within the set of passive states between 
canonical and 0-canonical distributions: they give rise to 
the most and least stable states, respectively. 

IV. APPLICATIONS 

Besides the question of activation [8, 18, 19], there 
are other scenarios related to work extraction where the 
MBPS can be useful. 

Weight. The notion of passivity has been recently used 
to define the work in fully quantized heat engines [20-22] . 
In such set-ups, every component of the engine, includ¬ 
ing the weight which serves as a receiver of the extracted 
work, is a quantum mechanical system. Then it becomes 
natural to interpret any change in Wmax, the work con¬ 
tent of the weight, as the work exchanged with the engine 
[20-22]. 

In order to relate the work stored in the battery with 
the standard notions of thermodynamics, such as the 
Carnot principle, it is useful to obtain bounds on (1) 
that do not depend on the whole spectrum of the state, 
but rather only on its energy and entropy. It is straight¬ 
forward to see that (1) satisfies, Wmax < tr((p — Tpi)H) 
where S{p) = The authors of [20-22] use this 

bound to obtain an upper bound on the Carnot effi¬ 
ciency in fully quantized set-ups, where the entropy gain 
of the battery is non-negligible. Now, using the MBPS, 
we can find a bound on Wmax in the other direction, i.e., 
Wmax > tr((p - crpiJ), where S{p) = 5'(ct*). This ex¬ 
pression provides a lower bound on the extractable work 
from a battery, which, again, only depends on the energy 
and the entropy of the state. 

Thermal bath as an ancillary system. Passive states 
can also be activated if one has an access to a ther¬ 
mal bath at some inverse temperature /3. Then the sec¬ 
ond law of thermodynamics places an upper bound on 
the extractable work through the free energy difference 
of the system, W < Fp[ap] — where Fp[p] = 

tv{Hp) — j3~^S{p) is the free energy, and oc e~^^ is 
the thermal state (see, e.g., [3, 23-26]). The inequality 
can be saturated if the thermal bath is big enough and 
it is capable of thermalizing the system [26]. 

iS) 

Let us now define a fictitious thermal state, , whose 
temperature is adjusted to satisfy S'(T^f^) = S'(crp). We 


can now express the free energy difference as 

Fp[ap]-Fp[Tf] = /^^US,Ep)+F{T^p^^)-F{rf). (9) 

As in the previous case, this expression is maximized for 
Op = tJ*, and thus a* provides a bound on the amount 
of work that can be extracted from a passive state using 
an external bath. 


V. SPECTRUM 

The amount of work Amax(5'j A*) that can be locked 
in (T* highly depends on the structure of H and its di¬ 
mension. As an extreme case, when the dimension d of 
the system is 2, all passive states are thermal and thus 
£^maxiS,E*) = 0. As the dimension increases, so does 
Amax(<S', Bp), with a rate defined by the structure of H. 
In this section we give some general considerations in the 
limit of d —>■ oo. These asymptotic results are then illus¬ 
trated by exactly solving some specific systems for finite 
dimensions. 

Subexponential growth of the density of states with en¬ 
ergy. Let us assume a dense spectrum bounded from 
above by Em (the ground state is taken to be non¬ 
degenerate and to have zero energy). Assume that the 
density of states (DOS) [27] scales polynomially with 
energy, g^^ = cE°‘, where c is some positive constant. 
The total number of states within [0, E] is then given 

by Ag = JgdE'g^, = Let us define as a 

state that is filled up to energy B, i.e., in 

(6). It satisfies = ■^J^dE'g^,E' = ^B, and 

= InAg. The MBPS is a combination of two 
such states, -I- (I — A)a;g^, with Bi, B 2 depend¬ 

ing on the specific case (entropy of the state, spectrum, 
etc). Numerical analysis provides evidence that Bi = 0 
and B 2 = Bni is always the optimal choice for 1. 

Therefore we focus on ctq = (1 — A)|0)(0| -|- Xuje^, where 
A is determined by the energy (or entropy) of tJo- 

The energy and entropy of ao are given by El^a^) = 
tr[croB] = A^^Bm and S'(cro) = B(A) -I- AlnA^j.^^ -|- 

O , where H{X) = —A In A — (1 — A) ln(l — A) is 

the binary entropy in natural units of information. From 
B((To) and S{ao), one can express the entropy as a func¬ 
tion of the energy, S{E). Notice that S{E) —>• 0 for 
in^rn. pj That is, if the norm of the Hamiltonian, 

Em, is big enough, then essentially the state has zero en¬ 
tropy while having a finite energy. This is in sharp con¬ 
trast with a thermal state, where if 5” —)■ 0 then B —)■ 0. 
This also implies that the energy hidden in a passive state 
^max('S', E*) can be arbitrarily large in the d —>■ 00 limit. 

Bathlike spectrum. Assume that the DOS now scales 
as = e^^ . Focusing again on ao, we obtain E{ao) = 
A (Em — b~^) O (e“^^™). From this expression we de¬ 
termine A(B) which, together with Ne = (e^^™ — l)/b, 
can be inserted into S'(cro) to calculate S{E). Taking, 
again, the limit E/Em —>■ 0, we find S —)■ bE. Here 
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FIG. 2. Aniax(5', i?*)/lnd versus 5/lnd (a) for an equally 
spaced Hamiltonian with d = 50,100,200,400; (b) for a 
collection of n non-interacting two level systems with n = 
10,50,100,200 (and d = 2"). As the dimension increases, in 
(a) so does the energy difference between the most energetic 
passive state and the thermal state, while in (b) the difference 
grows much slower due to the presence of large degeneracies. 
Insets: Amax(‘S', i?p)/Ind versus d or n, for fixed small value 
of S. While in (a) there is linear growth, in (b) the value 
grows only logarithmically. 

we observe that any amount of energy has an associated 
amount of entropy, even if the Hamiltonian is unbounded 
from above. Thus here, contrary to the previous case, the 
energy has a heatlike character. 

To illustrate this behaviour, in Fig. 2 we compute ex¬ 
actly Amax(5', E*) for a single system with equally spaced 

eigenvalues, H = for a collection of n 

non-interacting two-level systems [28]. It is clearly ob¬ 
served how the presence of high degeneracies hinders the 
growth of Amax(-S',i?p) in the latter case, while in the 
former the energy difference grows with system size. 

Finally, we note that an exponential scaling of the DOS 
with energy is a common assumption for the Hamiltonian 
of a bath. Indeed, a thermal bath is expected to be stable, 
has to thermalize any system that is brought in contact 
with it, and has to be a large system with macroscopic 
(extensive) energy and entropy. On the other hand, the 
stability [29], extensivity and the ability to thermalize 
[3, 11-14, 30] hold only for systems with short-range in¬ 
teractions, which, in turn, have exponential DOS. Our 
results can thus be seen to provide a new insight on the 
role of the DOS. As opposed to polynomially growing 
DOS, passive states with exponential growth appear to 
behave pretty much like standard thermal states. This 
is in the spirit of the equivalence of canonical and mi- 
crocanonical equilibria, that, again, holds only for sys¬ 
tems with short range interactions [3, 11-13]. It is worth 
adding that these type of spectra play an important role 


in fundamental questions such as thermalization [30] or 
the third law [31]. 


VI. CONCLUSIONS 

In this work we have characterized the family of passive 
states that maximize the energy of a system for a given 
entropy. We proved that they also solve the dual problem 
- they minimize the energy for a given entropy. There is 
thus a clear parallelism with thermal states, which pro¬ 
vide the reverse solution to such optimizations. These 
extremal properties make this class of states useful to 
obtain bounds in quantum thermodynamics of finite di¬ 
mensional systems. Indeed, we have shown that this class 
provides a lower bound on the amount of work that can 
be extracted from a thermally isolated quantum system; 
and it places upper bounds on the extractable work from 
a set of passive states. 

We have also discussed how energy and entropy are 
related for the MBPS depending on the spectrum of the 
Hamiltonian. Whereas in thermal states any amount of 
energy is associated with some gain in entropy, we have 
shown that this is no longer true for (8) if the spectral 
density of the Hamiltonian increases sub-exponentially. 
This demonstrates a clean cut between bath-like spec¬ 
tra (collections of systems interacting with short ranged 
forces) and other types of spectra (systems with long 
range interactions). 

Finally, the family of states found here can be used to 
lower bound the extractable work from a set of correlated 
states, complementing the results in [32]. We leave as 
a future work to further explore the implications that 
the MBPS have for the efficiency of fully quantized heat 
engines [20-22], for generalized notions of passive states 
[33, 34], and for other scenarios where thermodynamic 
processes are modelled by unitary operations [24, 35-40]. 
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